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1 Introduction 



One of the outstanding problems in the theory of two-dimensional quantum gravity 
is the effect of matter fields on the Hausdorff dimension. In models of discretized 
two-dimensional quantum gravity we define the grand canonical partition function 
for an ensemble of graphs (which for the moment we assume are triangulations) Q 
by 

^(/^) = E^"'""^G (1) 

Gee 

where |G| denotes the number of triangles in G, and wq the partition function of 
any matter fields in the theory on the graph G (for an introduction to this material 
see for example To define the Hausdorff dimension 0, Q we first define the 
geodesic distance dG{i,j) between two links i and j as the minimum number of 
triangles which must be traversed to get from the centre of one link to the centre of 
the other. Then we introduce the two-point function 

n{r,f,) = ^e-^I^I^G SidG{hj)-r). (2) 

We expect that H has the asymptotic behaviour [|], || 

n{r,fi)r~. e-™(^)^ m(^)r » 1, 

~ r^-'^s^ m(/i)"^ >> r >> 1, (3) 

where, as — fj,c, the mass gap vanishes as 

m{fi) ^ ifi- fi,Y^. (4) 

In general it is also convenient to consider a more general correlation function be- 
tween boundary loops of length li and I2', (0) is essentially the correlator for minimum 
length loops. Note that it follows from (^ that 

5^H(r,/i) = 5^e-'^l«l|G|'u;G 

r Gee 

~ (/x-/x.)-^- (5) 

where, in unitary theories, 7str is the string susceptibility exponent, and inserting 
the form (||) we conclude that 

^9(2 - ^g) = 7str (6) 

which is the Fisher scaling relation. At least in unitary theories the Hausdorff 
dimension dn is given by dni^g = 1 and has the geometrical meaning that in the 
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continuum limit the average volume is related to the geodesic size by {V) ~ W^" 

i- 

Analytic calculations of the scaling behaviour of the correlation functions 
were first done by means of the slicing decomposition introduced by Kawai et al 
and then somewhat later Watabiki ||^ introduced the peeling decomposition. For 
pure gravity (ie = 1, 7str = both peeling and slicing decompositions keep 

track of the geodesic distance and give the same results which tell us directly that 
the Hausdorff dimension of the ensemble is 4. 

When matter fields are introduced the situation becomes more complicated. The 
time scale, usually called the string time t, introduced in the decompositions that 
have been formulated is no longer by construction the geodesic distance, nor indeed 
are the time scales for different decompositions necessarily equivalent. However the 
above discussion of correlation functions can be repeated in terms of the string time 
t instead of the geodesic distance r leading to another pair of exponents, r] and z/, 
which are also expected to satisfy the Fisher scaling relation. For example in the 
c = —2 model the scaling with string time has been calculated completely by the 
peeling decomposition with the result that = \ which would imply that dn = 2 
if the string time and geodesic distance are proportional. In fact high precision 
numerical calculations find dn = 3.58 ± 0.04 in agreement with the formula 



dn = 2 ==■ (7) 

V25 - c + vl - c 

derived using scaling arguments for diffusion in Liouville theory. For unitary matter 
complete calculations have not been made but it seems that z/ = |7lsti./2 0; the 
implied value of d^ is in contradiction with the results of numerical simulations 



which suggest that dn is close to 4 |jTO[ but are not in particularly good agreement 
with either. It seems certain that when matter is present the string time and the 
geodesic distance have different scaling dimensions but the relation between them 
is unknown. 

In this paper we will be concerned with the matrix models with a critical point 
corresponding to the (p, q) = (2, 2K — 1), K = 2,3,4, . . . multi-critical models cou- 
pled to quantum gravity |Tl|, [1^. For K > 2 the ensemble of graphs Q now allows 
polygons with {4, . . . 2K} sides and there are K independent coupling constants, 
with polygons of order {6, 10, 14, . . .} having negative weights which makes the mod- 
els non-unitary. The partition function is again defined by 

Z(A.) = 5^e-^I^W (8) 
Gee 

where now G denotes the number of polygons, fi is the coupling constant conjugate to 
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the number of polygons, and wq depends on the remaining K — 1 couphng constants 
of the theory. K = 2 corresponds to pure gravity; the second couphng constant 
which appears in this matrix model is conjugate to the length of the boundary of the 
graph. Since there are now in general many couplings, there are many correlation 
functions which are second derivatives with respect to the couplings so there are 
many susceptibilites. We define the susceptibility 

~ (/i-/ic)-^ (9) 

where the exponent 7 is known to take the value —K~^ at the multicritical point. 
The string time, tQ{i, j), separating two links is now defined as the minimum number 
of polygons which must be traversed to get from the centre of one link to the centre 
of the other and the two-point function is defined as 

= ^e-'^l^ly^G Yl ^(^g(^,j) -i)- (10) 

Gee i,j6G 

We expect (and shall confirm) that, if all the couplings except fi are set to their 
values at the multi-critical point, Ti. has the asymptotic behaviour 

nit, fx) ~ e-'"(^)*, m{fx)t » 1, 

~ t^-^ m{ixy^»t»l, (11) 

where, as yU ^ /ic, the mass gap vanishes as 

m(/i) ~ (yU - ficY- (12) 

and that the Fisher scaling relation (which may be obtained by similar manipulations 
as before) 

iy{2-v)=7 (13) 

is satisfied. One could hold that t is still the geodesic distance but this is slightly 
problematic for large K; it implies that all sides of a given polygon, no matter how 
many sides it has, are separated from one another by geodesic distance 1 and we shall 
argue in section 6 that the continuum limit of t is not a sensible continuum geodesic 
distance. These models have been analyzed using the slicing decomposition in the 
same way as pure gravity |T^ and they have also been considered using the peeling 
decomposition in the scaling limit but the discretized equations have not been solved 
completely ||^, |T^. In this paper we will examine their peeling decomposition in 
detail and explain how to solve completely the non-trivial differential equations that 
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Figure 1: The Schwinger-Dyson equation for the potential (0). 

arise. We have two motivations for this; to check whether the results are indeed the 
same as for slicing, and the intrinsic interest of the method of solution. 

This paper is structured as follows. In section 2 we briefly describe the standard 
peeling calculation for pure gravity and then derive the evolution equation for the 
multi-critical models. In section 3 we consider the K = 2 case and show that it 
always gives the standard pure gravity results. Then in section 4 we show in detail 
how to calculate the i] exponent for K = 4 and describe how the calculation extends 
to all higher even K. In section 5 we explain how to calculate u for all even K and 
in section 6 we give our conclusions. 



2 The Peeling Decomposition and Evolution Equa- 
tions 

We start by reviewing the calculation in |^ for the simplest pure gravity model 
which has matrix model potential 

f/(0) = - Ig^'- (14) 

The matrix model generates the dual graphs to the triangulations Q (see equation 
(1)) with g = and wq = 1- The Schwinger Dyson equation for connected Green's 
functions is obtained by marking one external line and pulling it out to expose the 
vertex to which it is attached |jl5|, see fig.(p. In the peeling decomposition we assign 



a time variable to this process; a single iteration advances t by an amount 1/n so 
that we obtain 

n 

n—m+1 it). (15) 

m=l 

We are interested in the loop-loop correlation function; suppose for the moment that 
at t = the entry loop is a one-loop and form the quantity 
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Figure 2: The two-loop function in spherical topology. 



which is the amplitude for an exit n-loop at time t given an entry 1-loop at time 
t = 0. Differentiating (|T5|) we obtain 

n 

Gn{t + 1/n) = gGn+i(t) + 2gY, Am(t)G 

n—m+l it). (17) 

m=l 

If we restrict to spherical topology then Am{t) may be replaced by the disk amplitude 
for m legs, Am, because the branch can never rejoin the main tube (see fig.(^). The 
next step is to approximate the time by a continuous variable to obtain the evolution 
equation 

+ 1^^^ = gG^+,{t) + 2gf2 A^G„_„+i(t) (18) 
n at ^-^ 

m=l 



with the initial condition that 



- '-urn ^ <^«' 



Note that by differentiating (^8|) we can show iteratively that all the derivatives of 
Gn{t) at t = are finite. Defining the generating function 

oo 

G{t,x) = J2^''Gn{t) (20) 

n=l 



9G arn^^ 



equation ( [T8| ) becomes 

dt ~ dx 
where 

F{x) = 2gA{x) +g-x, (22) 
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and the form of the disk amphtude 

oo 

A{x) = Y,^nX^ (23) 

n=l 

is known Note that the function F{x) contains just the universal scahng part of 
the amphtude A{x)] it is given by 

F(x) = (/-x)(l-8/(7-4(7x)^ (24) 

where / is the root of the cubic 

(l-8/^?)f -/ = (25) 

which is positive and vanishes as 5^ — > 0. It is then straightforward to solve the 
evolution equation which gives 

G{t,x) = ^U{t + J{x)) (26) 
r [x) 



where 

dJ 1 



dx F{x) 

and the function U{y) is fixed by the initial conditions to satisfy 

1 



(27) 



F{x_ 

This leads to the scaling behaviour at large y 

cosh(/i - ^c)^y 
sinh^(/i - iJc)iy 



U{J{x)). (2J 



^(,) ^ (;.-/..) -^-(7 (29) 



where we have suppressed various constant factors; so we deduce = \ and rj = A 
in agreement with the Fisher scaling relation 0. 

Now we turn to the multi-critical models for which we will use the notation of 



reference |TH]. The potential for the i^-th multi-critical model is given by 

K 
p=l 

where the couplings at the multicritical point are 

- _ U . (31) 
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Figure 3: The Schwinger- Dyson equation for disconnected graphs of the multicritical 
potential (|30|) . 

This time we will deal with disconnected graphs. The evolution equation follows 
from the Schwinger-Dyson equation, shown in fig.^, just as in the pure gravity case; 
we obtain the evolution equation 

n-2 K 

A^t + 1/n) = J2 An-j-2{t)Aj{t) -2J2 kg^An+2k-2{t). (32) 
Proceeding as before, but assuming that the entry loop is an m-loop and defining 

we find 

Gn,m{t) + — ^ ^ = 2 Gj^rn{t)An-j-2 — '^kg^ Gn+2k-2,m{t) (34) 

j=0 k=l 

with the initial condition that 

(0) = Sn,m. (35) 



Again, note that by differentiating (p^ ) we can show iteratively that all the deriva- 
tives of Gn,mit) at t = are finite. 
Defining the generating function 

oo 

GUx,t) = J2^''Gn,m{t) (36) 



n=l 



we obtain the partial differential equation 



K 



dGm <9 J ( o \ 1 2A;c/f 



K 2k-3 

2%f E GUtm - 2 - j)x^-''^+\ (37) 
k=i j=i 
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At this stage we should make several remarks. Firstly that one reason for dealing 
with the disconnected graphs is technical convenience; the required disk amplitudes 



are known and take a simple form , and the structure of the evolution equations is 
similar to the 0'^ case. However it is also very clear that because the graph ensemble 
is disconnected there is no direct correspondence between the string time t and a 
geodesic distance - the latter can only be sensibly defined on connected graphs. 
These equations are more difficult to solve than the pure gravity example reviewed 
above because of the presence of the a priori unknown functions Gj,m{t) which have 
to be determined by the self-consistency and analyticity properties of the solutions. 
We will work our way through the problem in a number of steps. First we will study 
the solution in the K = 2 case and show that it gives the standard pure gravity 
results. 



3 Universality in the K = 2 case 

It is simpler to work at the multi-critical points initially and to compute the rj 
exponent directly at the critical point; we leave the v exponent for section 5. For 
the K = 2 multi-critical model we have 

BC f) 

^ = X— {Gmi2x'Aix) - 1 - 2gl - Aglx-') - Aglx-'G^,^{t) (38) 
where gl = 1/2, = —1/12 and A{x) is given by 

Thus Cm satisfies 



dx-^G^ d /G^fl- 4x2)1 \ 1 



dt dx \ 3a;2 / 3^^ 
with the initial condition 

GUO,x) = x"'. (41) 

The potential (^) is an even function of the fields and therefore Gn,m{t), which 
is the amplitude for an entrance m-loop and exit n-loop, can only be non-zero of n 
and m are both odd or both even. Thus if m is odd there is an unknown function 
on the r.h.s. of ( ^OD whereas if m is even there is no such problem; we will consider 
the even and odd cases separately. 
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3.1 Even m 

When m is even Gi^m{t) = and we can solve (|iO|) immediately to obtain 

GUt,x) = ^^U{t + J{x)) (42) 

where 

fl -4x2)1 



F{x) 



3x 
dJ 1 



dx Fix) 



J{x) = ^(il-AxY^ -1), (43) 



4 

and the function U{y) is determined by the initial condition (^ 



Thus we obtain 



for large y and hence 



X- = ^fiAx)). (44) 

Uiy) ~ 4 (45) 

G„,„(t) ~ ^ (46) 
for even n and m so that = 4 as expected for all these amplitudes. 

3.2 Odd m 

When m is odd G'i,m(^) 7^ and the presence of the unknown function on the r.h.s. 
of (|40|) complicates matters; however this is more typical of the general multi-critical 
models than the even m case and so we shall study it in some detail. Although the 
differential equation can of course still be solved in the t domain it is more convenient 
to work with the Laplace transformed correlation functions 



poo 

Gn,m{s) = / Gn,m(^)e ^ dt 

Jo 



CO 



GUs,x) = - I GUt,x)e-''dt. (47) 
X Jo 

Taking the Laplace transform of (^O]) we obtain the equation 

sGm{s,x)-x^-' = ^ {G^{s,x)F{x)) + ^Gi,^(s). (48) 
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Integrating this differential equation we find 

^ ' F{x) V 3x io V (i-v)i;; ^ ^ 



The function Gm{s,x) has a power series expansion in x 

= ^x"-iG„,„(s) (50) 

n 

which we expect from (|^) and (|43|) to have finite radius of convergence |; within the 
radius of convergence the coefficients are the Laplace transforms of the correlation 
functions. Unless the Gn,m{t) grow faster than exponentially at large t, something 
we do not expect, their Laplace transforms Gn,m{s) will have an asymptotic series 
representation at large s. This series can be obtained by successive integration by 
parts of the definition (^); as we observed in section 2, Gn,m{t) and all its (finite 
order) derivatives are finite at t = so we obtain the formal series 



Gr, 



"o = <5„,m. (51) 



s ^ — ' s 

k=0 



(Of course the ak depend on n and m but we will always suppress such dependence 
for clarity.) We will now show that imposing the condition that the coefficients 
in the x expansion of Gm{s,x) behave like (^TJ) at large s is sufficient to fix the 
unknown function (ji,m(s). 

It is more convenient to impose the condition on the integral of Gm.{s,x), 



GUs, y)dy = e''^-^ / dye-^'^y^ y^~^ - ''""l , . (52) 

sm S ./n \ (1 - 4?/2)- 

Using (13), (0) and ( |52D the consistency condition is that 



,SJ,.J I ^y^-sj,y) I ym-^ ^ (53) 

(i-4?/^)2y s 

That is to say the small x, large s, expansion of Qm contains no terms 0(1) or higher 
in s. Substituting (|51| ) and expanding the exponential in its Taylor series (which is 
of course allowed for all values of the argument) we get 
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For x < I (^4|) clearly has a power series expansion in x so the requirement that 
there are no 0(1) terms in s becomes 

l-X f ^ oo \ 

C,= \ dy[ y-~' - ( J(x) - J(2/)r J = 0. (55) 

First we show that this constraint alone implies that all higher powers of s vanish 
as well; for 0(s) we get 



C, = [ iy ( - J(y)) [y-^ - ^^-i-^ g (J(.) - J(y)f ^ j = 0. 

(56) 

and observe that 

^ = = (57) 

ax ax 

provided x < |. It follows that Ci{x) is a constant; but Ci(0) = therefore 

Ci{x) = 0. (58) 

Any pair Ck{x) and Ck+i{x), with k > are related in the same way and so it is 
straightforward to proceed inductively to show that all Ck>o{x) are zero. 
Taking the Laplace transform of (^) with respect to J we find that 

-e-'-Uj = Gi,Us) ^Tl^dJ (59) 



m ' 7q ?>x{J) 

where x{J) is obtained by inverting (^). Thus Gi,m(s) is determined; making the 
change of variables x = itanh0 and integrating by parts we get 



tanh™-^ e-i^^°^^'^sech2' 
2'"-is/~cosh(/.e-t^ 



It is clear that G'i,m(s) is positive for all real, positive s and straightforward to 
check that its second derivative with respect to s diverges logarithmically as s — > 0; 
it follows that 

GiAt) ~ ^ (61) 



at large t. The remaining integrals in (^) yield functions which are analytic in s in 
the neighbourhood of the origin provided x < Xc- Thus we can conclude that 

G„,m(t) ~ ^ (62) 

for all odd n and m. So the result rj = 4 is universal for all amplitudes Gn,m{t) in 
agreement with every other calculation for pure gravity. 
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In fact can be fixed more simply by examining (^). As x | Xc, J{x) 

diverges and lience Gm{s, x) grows faster than exponentially in s, which is impossible, 
unless 

dxe-''^^^ ( x^-' - \ = (63) 



[l-Ax^y. 

which is the same condition as (pUD- We have explored the more indirect route 
because this will help in the multi-critical case. 



4 The Multi-critical Models 



The multi-critical evolution equation (|37|) may be written 
dx^^Gmix-it) d 



dt 



dx 



{x-'Gm{x,t)F{x)} 

K 2fc-3 

- 2%f ^ GUmk - 2 - j)x=-''^' (64) 
k=i j=i 



where, taking A(x) from [0], we have 

F{x) = ^ ■ ^ 

\ k= 

= -2Kg^x'-'''{l-AxY--^. 



(65) 



The solution to (p^ ) takes the form 
1 



F{x 



-Uo{t + J{x)) 



E2A:^?f X;^^^^Y^ fdrG,,rn{r)U,,{t-T + J{x)) (66) 
k=i j=i ^ ' 



where, as before, 



J{x) 



dy 
F{y) 



K-2 



^/^f(l-4x 



K~l 



(67) 



r=0 



and the coefficients are easy to compute. The function Uo{y) is fixed by the 
initial conditions 



x 



m—l 



1 



Fix 



-Uo{J{x)) 



and the functions Ujkiy) by requiring that (|^) is a solution of (16^), 



X- 



j-2k+l 



Fix 



-U,,{J{x)). 



(6J 



(69) 
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Because of the unknown functions Gj^mit), (|66D is not of course a complete solution. 
Note that J{x) is always an even function which diverg I — i> ^ and which is 

positive for even K and negative for odd K. To determine the critical behaviour 
from the properties at large t we need to know the behaviour of the functions Uo{t) 
and Ujkir) for large positive argument. For even K we see by considering (|68|) and 
iW^ as la; I 



I that 



2K-1 



Uq[t) ~ r 2^f-3, 
f^.fcM ™ (70) 



However for odd K only the large negative argument behaviour is determined. This 
phenomenon always occurs in calculations for the multi-critical models. Extrapolat- 
ing the solution to positive time leads to a singularity at finite time; it seems to us 
quite likely that this is an artefact of the truncation of the original finite difference 
equation (|3^) into a first order differential equation p^ ) and that the solution may 
be stabilised by higher derivative terms. From now on we will concentrate on the 
even K models and start by studying the K = A case in detail. 

4.1 K = A, even m 

As for the K = 2 model we will consider the cases of even and odd m separately. 
For even m we have 

dx-'GU x,t) d ,^ , ^ / 4 

dt 

with 



^{x-^G^(x,t)F(x)}+G2,™(t) - + ^G4,„.(t) 

(71) 

= ^ o. V • (72) 



Taking the Laplace transform gives 

sGUs,x) - x"'"^ = ^ {Gmis,x)F{x)) +G2,Us) - ^) + ^^^'-("^ 

(73) 

where the transformed correlation functions are defined as in (^); integrating ( |75D 
we obtain 



+ ^^^4(^(y'Ge':^(^) +y'Gj„(.) +i/G2,™(s)))| (74) 



13 



where we have introduced the combinations 

gI^{s) = 70G2,™(s)-14G4,„(s) + G6,™(s). (75) 

Note the appearance of Ge.mls) in (0); this happens because F{x) is singular at 
X = which makes the evaluation of the limits of integration slightly non-trivial. 
As before it is convenient to deal with 

Gm{s,y)dy= h— (76) 

sm s 

where now 



35i/5F(y) 

and we deduce that 



{y%^{s) + y'Gl^is) + yG,,Us)) } (77) 

^7™ ~ - (78) 
s 

in order to fulfil the constraint that Gm{s,x) has the correct small x, large s, ex- 
pansion. In fact we can immediately determine GQ^m{s) as explained at the end of 
section 2 by requiring that Gm{s, x) does not grow faster than any exponential of s 
as a; t I which implies that 

9m{s, \) = (79) 

and so GQ^m{s) is related to G2,m{s) and G4,m(s) which still have to be determined. 
Defining 

Ge^mis) = -^-^^ 

k=0 
A:=0 

G^A^) = -Y.\ (80) 

k=0 

(of course the a, [3 and 7 coefficients depend on m but we have suppressed this to 
avoid clutter) and substituting into (^) we obtain the 0(s°) constraint 



0. 

(811 



14 



Identical manipulations to the K = 2 case show that Cq = ensures Ck>o = 
also. There are in fact three separate constraints hidden in (^Tj) which are sufficient 
to determine the a^, Pk and 7^ coefficients. To see this we proceed by writing 
everything in terms of J. J{x) is an even function with a small x expansion 

J{x) = f x« + + . . . (82) 

o 4 

so it follows that, by reverting the series, 

^l:^ = (jl)W2^^^(j|)fc (83) 

k=0 

where we take the real positive cube root of J (remember that m is even and again 
we suppress the m dependence on a^). Similarly we have 

dJ' 



/ dj'{j-j')'j'--^j2^i{j'-^y 

Vfc,jfc+(m)/3 r(fc + l)r(^) 



and 



^5) 



1=0 - I 3 

We will not need explicit expressions for the coefficients ai, hi and q. Substituting 



and (|85D into Co we obtain 

00 00 ■p(l+2\ 

+EEA^^^^^""" ^(/'il,' V (86) 



fc=0 «=0 V ' 3 



N-\ 



Letting N be an integer, ([86| ) yields three conditions corresponding to terms 0( J 
0(7^"*"^) and 0{J^^i) respectively. The corresponding coefficients are easily ex- 
tracted from (^) and we find 

N N 
- a3^-„/2r(iV + 1) + ^ ^M-LhL-lT{L) + /57V-LC3L-2r(L) = 0, (87) 

L=l L=l 

15 



4 ^ 1 ^ 1 
a^N+i-^/2T{N + 3) + 5^ -iN-LhLT{L + 3) + 5^ /9jv-LC3L_ir(L + 3) = 0, (88) 

L=0 L=l 

5 ^ 2 ^ 2 

r(L + -) + 5^/3^„iC3Lr(L + -) = 0. (89) 

L=0 L=0 



Now suppose that 70 . . . 1n-i and /5o . . . are known; from (|87|) we can obtain 



ao . . . oat-i- Then (|88|) contains only these known coefficients together with 77V which 
is thus determined. Now ( PPj ) determines [3n- Knowing 70 . . . 7Ar-i and (3q . . . Pn-i 
we now determine from (0). Proceeding in this way all the coefficients, and thus 
the right hand sides of (|80D , can be obtained iteratively. To analyse the asymptotic 
behaviour of the correlation functions we need to know the behaviour of G2,m{s) etc 
at small s. Resumming the relations (|87|), (|88|), and (|^), and using the integral 
representation 

^^""^ ^^d^C-'e-^', (90) 



we get 



^ "'0 Jo ]^^^ 

+ Kmis) rf^e-^^ 5^ C3^-2e'^-\ (91) 

"^0 7V=1 

/OO ^ POO ^ 

poo 

+ Glmis) ^^ee-«^5^C3^-ie'^-^ (92) 

■^0 7V=1 



= - / f^ee-«^$^a3iv+2-W2e^^^+G2,„.(s) / rf^e-^^ J] We^"^ 

/•OO ^ 

+ Gi;^(^) / ci^e-^^ 5^ C3;ve^-i (93) 

These three equations determine the unknown functions that we need. 
Now consider the equation 

J{x) = e (94) 

where ^ is real and positive; it can be rewritten as the quintic equation for u = x"^ 

24 

(1 + yO'(l - 4u)' - {30u^ - lOu + If = 0. (95) 
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When = it reduces to 



- 4m^(10 - 95m + 25Qu') = (96) 

which has a complex conjugate pair of roots 

95 ± i9VT5 , , 

and three roots vanishing when = 0; any integer power of these roots, which we 
denote by Mi 2 3, has a series expansion 



n=0 

00 

n=0 

00 

{usiOr = {u:'i^rY.DMn{u:\^T. (98) 



n=0 

1 



where is a complex cube root of unity and .^3 is the real cube root of ^. 

Now take linear combinations of (pT])- ( pSf ) with the coefficients being different 
powers of uj to obtain 

with i = 1,2, 3. 

As ^ cxD all the roots of (|95D converge on u = ^ like 

u = ](l-a^ + ..) (100) 



4V 

where a is a fifth root of unity and c is a constant. The fiow of the roots in the 
complex plane as C, varies is shown in fig.§. We see that ui fiows into the C, = 00 
point as ( |100| ) with a = 1; and that the complex conjugate pair ^2,3 fiow into the 
pair with a = e^*^. Note that there are no degenerate roots of (0) in the interval 
^ = (0, 00) (this follows from the fact that F is regular in this interval). 
Now define the function J' by 

J(u = x^) = J(x) (101) 

and make the change of variable 

e = J{U^) (102) 
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Figure 4: The flow of tlie roots of ( P^D in the complex u plane; the arrows denote 
the direction of increasing ^. 

in (p9|). The constraint equations then take their final form 

where the contours Pi in the complex w-plane are shown in fig.(^). Note that the 
i = 1 constraint is the same as ( [79D after an integration by parts and that these 
equations are guaranteed to have a unique solution by the argument immediately 
following ©-(ID. 

To find the asymptotic large t dependence of the correlation functions it suffices 
to find the leading non-trivial s-dependence of their Laplace tranforms at small s; 
for this purpose we need the integrals in (|103|) up to and including O(s^). Replacing 
by its explicit form we need 

If = ^-f duu\l-Au)-''/^J^{u)exp{~s{A{l~Au)-^+B{l-Au)-l + C{l-Auy"^)) 
2 Jp. 

(104) 

where J^{u) = jn^^ u^™}, A = 7/64, B = -35/96 and C = 35/64 and we 
have dropped a factor exp(7s/24) which cancels in ( |103| ). Now make the change of 
variables 

z = s{l - Au)-^ . (105) 
Then all the integrals in ( p.02| ) are linear combinations of 

= J dzz-^^^~^^/^ exp{-Az) exp (^-{Bshr> + Cst^i)) (106) 
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Figure 5: Integration contours for ( p.06|) in the complex z plane. All contours 
start aX z = s and go out to infinity; the cut is for the fractional powers of z. The 
contour Pg is simply the complex conjugate of Pg- 

where k = 1,2,.... The contours P/ in the complex 2;-plane are shown in fig.^. 
The second exponential in ( |10(j| ) is then expanded in its Taylor series, and the P2 3 
contours collapsed onto the branch cut as shown in fig.^; all the integals then become 

A A 




Figure 6: a) Distortion of the integration contour P2. If the integrand has an 
integrable singularity at the origin then the circular part of the new contour can 
also be collapsed onto the branch cut as shown in b). 



elementary. The equations ( |103| ) can then be solved to give 



GA,m(s) 



+ 2) + ^:^f7fy(^ + 6m + 8)(1 + 2 cos -) ( - j + 0(.^), 



35 

26+m 

35 



175 rf-) 

(3m-2) + ^^(m2 + 6m + 8)(l + 2cos^) (^) +0(.t). 



V 64 



26^ 



-(15m2-50m + 64) 



875 r(t 



27+m p( 



|((m^ + 6m + 8)(l + 2cos- 



47r, /7s\ 5 



5 ' V64y 



+ 0(si 



(107) 



These expressions have the expected properties; they are positive and decreasing 
functions of s for small positive s. It follows from ( |107|) that 6*2,™ (^), G4:,m{t) and 



19 



Ge^mit) all have the same asymptotic behaviour at large t namely 

G,,Ut) ~ t"^ (108) 

for j = 2,4,6 and hence the exponent rj = 12/5. To determine the behaviour of 
the higher correlation functions it is simplest to return to (|7^. All the integrals 
yield functions of s which are analytic in some neighbourhood of the origin provided 
a; < 1/2 and so we can conclude that every coefficient of the x expansion of Gm.{x, t) 
behaves the same way and that ( |108| ) is valid for all (even) j and m. 

4.2 K = 4, odd m 

The method is very similar to the even m case. After Laplace transforming the 
evolution equation and solving for Gm{s,x) we obtain 







+ 35/ir(^) {y'Glm{s) + y'Gl^{s) + Gi.^(.))) I (109) 
The consistency conditions can be cast in the form 

= /d„^e---<«.|„l'»_l(£l^ + S=lW+^H ,U0) 

Jpi du I 35 1-^2 U2 U2 I \ 

This time we need the integrals 

If = ^[ rfn4(l-4^/)-'/'{^/-^n-^l,M^(™+l)}e-^^(") (111) 



JP, U2 

To the required order in s these can be calculated by expanding the factor about 
M = |; then term by term the integrals are just the TCk- The resulting expressions 
for the singular terms in s are simple but the constants appear as infinite sums; 
fortunately the constants are only needed for the sub-leading s dependence so this 
method suffices (the reason for this is explained in the next sub-section). The final 
expressions for the correlation functions are very similar to those for even m 

35 

G,,mis) = ^^(m^ + 4m + 3) 

175 r(^)^ „ „ Att, f7sY^ ^, 4, 
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35 



175 r(^)^ , „ An, f7sY ^, 4- 

+ 2^fM^"^ +23m + 15)(l + 2cos-)(^-j + 0(.^ 



.5 

35 

G^A^) = T^zT— (15m2-20m + 29) 

375 Y(-) An /7s\ ^ 4 

+ ^^("^' + 9m2 + 23m+15)(l + 2cos— ) (^-j +0(si). 

(112) 

Thus we can conclude that the exponent rj always takes the value 12/5 in the K = A 
model. 

4.3 General even K 

The method is very similar to the K = A case. After Laplace transforming the 
evolution equation and solving for Gm{s, x) we obtain 

2K-2 



— If —c 

- 5JgfeE'^i''-'S^».(-'))} (113) 



where 



G2K-p-l,mis) = J2 Y^G2k-p-l,mi^). (114) 
k=2 



Note that for given m (either odd or even) half of the Gp^{s) in ( |113|) are auto- 
matically zero so there are K — 1 undetermined functions. Proceeding as before we 
next examine the roots of the equation 

j(x) = e (115) 

which can be rewritten as a degree 2K — 3 polynomial equation for u = x"^ taking 
the generic form 

(1 + aO'(l - ^uf''-' -{l + hu + ... + hK^2u'''^f = (116) 

where a, 6i, . . .bK-2 are constants. There are K — 1 roots which vanish at ,^ = 
and K — 2 roots which do not; the K — 1 vanishing roots correspond to the K — 1 
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functions which have to be determined. As ^ — oo all the roots converge on m = 1/4 
hke 

u = Ul-a^ + ...] (117) 

where c is a constant and a is a {2K — 3)rd root of unity. The roots which vanish 
at ^ = flow into 

f Ann \ K K , , 



From now on we will use the integer n to label the roots. The consistency conditions 
become 

» = / + 2A-9f E Jll^l (119) 

''P" I p=l J 

where the K — 1 paths P„ follow in the complex plane the roots which vanish at 
the origin. Note that in the case of i^' = 4 the conditions for even m, ( |103| ), and for 
odd m, (|110|) , can be merged into the form ( |119| ); for any fixed m we have K — 1 
equations for K — 1 unknowns. 



Restricting ourselves to even m all the integrals in (|119|) can be written as linear 
combinations of 

^P = J t^^i(l-4Mf-^+5exp -s^/if(l-4M)'^-^+ij . (120) 

Note that the factor exp{sh^_^) cancels out in (|119|) so from now on we drop it. 
After the substitution 

z = s(l - 4M)--^+i (121) 

we obtain 

2(2ir - 3)/Cp = s"^+^ y dzz'^ exp (^-h^z -Y^^h^s^z^'^^ (122) 

where the contours encircle the origin n times before heading off to real positive 
infinity, see fig. (7a). As we will demonstrate shortly all these integrals are needed 
up to and including the first positive power of s and so they can be separated into 
three classes; 

1. p < i^' — 1; the integrals are divergent as s but the singularities of the 
integrand are integrable at 2 = so the contour can be deformed as shown in 
fig. (7b). The integral from the starting point to the first blob is given by 

2/r — s 

2(2K - 3)/Cf " = ^^_^^^3 + Ois). (123) 
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b) 



Figure 7: a) A contour P' with n = 3 and b) the deformed version of it. The branch 
cut on the positive real axis is suppressed for clarity. 

The portion of the contour between the two blobs gives no contribution and 
the leg stretching out to infinity gives 



-div i27rn\ — IH — 



2K-3 



2(2X-3);C7 = {se 

^ / —i27Tn \ 2K-3 



dRR-^^ exp \-h^R - ^ h^R 



(124) 



To obtain the integral up to the desired order the second exponential factor 
can be Taylor expanded and then integrated term by term. 

2. p = K — 1] after an integration by parts the above construction can be used 
and we get 

2{2K-^)1Ck-i = (2ir-3) (^1 - (h^se-^'n^^ {^^)) ^ ' ' ' ' ^^^^^ 

3. p > K — 1] the leading non- analytic term is of higher order than we need to 
consider so 



2{2K-3)}Cp = — — ^TT^^. (126) 



2K-3 
2p-2K + 3' 

Now observe that in all these integrals s always appears with a factor e"*^'^" and 
that there are no other phase factors. Thus all of the constraint equations (|119|) can 
be obtained from the = case by making the replacement s —>■ se"*^'^". They can 
thus be written in the form 



K-l 

p=i 
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where 




TJw) =w''\y' Lnw^^ + Lkw + ...]. (128) 



We have truncated the expansion of J-'p{w) in anticipation of the following. The 
equations (|127| ) can now be written in the matrix form 



where 



nDf + s^fo^K-ii^ + 9LDfG + s^fG = (129) 



= exp(i27r(-ir/2+p)(-l + 2g/(2ir-3))) 



=pq 



fp — fo,p-l + fo,Ks5p^i, 

ujp = exp{i2n{-K/2+p)/{2K ~3)) , 

f = fp,q-l + fq,KS6p^i, 
=pq 

f = exp{i27i{-K/2+p)/{2K-3))fg,K-2 

=pq 

= ^2p,m- (130) 

Now D is clearly non-singular and it is straightforward to check that Q is non- 
singular. There does not seem to be any simple way of writing the elements of / for 
general K but we expect it too is non-singular. Then the leading order solution for 
G is a constant vector; furthermore the next term in the solution is 0{s'^^^'^^~^^). 
We have already shown by explicit solution that this is indeed what happens for 
K = 4] using Maple we have also checked it for = 6, 8, 10. It follows that 

V=— 131 

' 2K-3 ^ ^ 

Q 

for the correlation functions G^^^p = 2,4, ...2(i^ — 1). The relationship (|1 14|) 

Q 

between G^ ^^ and Gp_ Yn is non-singular so this conclusion applies also to Gp^rniP — 
2, 4, . . . 2{K — 1). Finally we note as usual that the integrals in (|113|) are analytic 
functions of s in the neighbourhood of the origin provided x < 1/2 so the conclusion 
extends to all the correlation functions. 



5 The z/ exponent 

To find u we need to study the scaling behaviour as the multi-critical point is 
approached and to do this consider the modified couplings 

9^ = 9^ 

^ = {1-A)g^, p>2. (132) 
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In the graphical expansion the power of 1 — A is the number of vertices in the graph 
so A is related to /x in (8) by 



1 — A = exp(/i — fic)- (133) 

By definition the multi-critical point is attained as A —>■ 0. The disk amplitude 
A{x) can still be calculated by exploiting the connection with topological gravity 
||T6|| ; the topological potential is 

V{z) = A~ KzA + z + {1- A){1- z)^ (134) 

and the disk amplitude 

A{A,x) = -l - 4^x')-i(l - V\z)) log ~ ^^^^ (135) 

Jc ^ 

where the contour encircles the branch cut of the logarithm. The branch points are 
at 2; = and 



2 = 1^1- tKA""" (136) 
where ex is a constant. Collapsing the contour onto the cut gives 

A{A,x) = K r'\y{l - Ayx\^^ (A + (1 - A)(l - yY'-') (137) 
^0 

Integrating we find that 

F{x) = x(2x^A{A,x)-l-J2^''^^ 



k=l 



-xKA{l - z^x^y^ - x(l - A) ^ ^(exAi/^)^-'=(l 

(138) 



which reproduces (|65D when A = 0. 

At small X, even for finite A we still get the leading behaviour 

F{x) ~ x^-^"" (139) 

which implies that 

Thus the considerations described in detail for K = 4 in section 4 up to equation 
(p3) go through as before - the only difference is that all the various coefficients 
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are now functions of A. However the contours Pi are modified; when we make the 
change of variable ( p.02| ) 

e = J{ui) (141) 

the ^ = cxo endpoint of the contour occurs at the value of Ui where J diverges. The 
points where J diverges are of course determined by the zeros of F. Defining w 
through 



X 



w 



(142) 



we get using (|138|) 

(cA-A^/^)^ r w-hiz^ + eK^^'^^w) 



■dw. (143) 



2 Joo KA + 2(1 - A)(eA'AVA')i^-i ^f^^ fc^f ^^-1 
The term KA in the denominator is sub-leading; its only effect is to shift the zeros 
by an amount 0(A^/^) and we discard it to obtain 



J{w) 



■ dw. 



(144) 



4(1 - A) J2ti kglf^'-' 

Note that the denominator has precisely K — 1 simple zeroes; one is real and positive 
and the others come in complex conjugate pairs. Each zero is the end-point of one of 
the K — 1 contours, see fig.(^. Since we will only need the leading scaling behaviour 

A 

P, 




Figure 8: The contours for K = 4 at finite A. 
we can also approximate {zc + exA^^^w) by Zc in ( |144| ) so that 

(6;,AV^)i-^ r 



A{l-A){-Kg§)J^ U,-^\w-Wk 



^'^^"^""y- YrJ -^dw (145) 



4(l-A)(-i^<7f)f^ V. w-w. 



k=l 
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where Wk is a zero of the denominator of ( |144]) and Rk the corresponding residue. 

Now suppose that w lies on the contour Pi which terminates at Wi; when w is 
close to the end-point we get 



-K 



Rl , ^/w + ^/Wl 

log ^ + Qi 



(146) 



where Qi is the accumulated finite contribution from all the other poles in 
Note that the branch cut for the square roots is taken down the negative real axis 
so that Rey/Wi > for all I. Inverting ( [146| ) we find that 

_ / r,exp(A'-»/'=''-g,J) + l V 
'"-""(^ r,exp(A--»/^«-g,J)-l j 

where Si and T/ are (complex) constants of 0(1). It turns out that the quantity 
Ri/V^i is never pure imaginary and thus ReSi ^ 0; except for the single real root, 
Si is not purely real either and hence the correlation functions have some oscillatory 
behaviour which is typical of a non-unitary theory. The formula ( p.471 ) is valid for 
J ^ oo] expanding we get 



w = Wi 



oo 

l + ^V'i'^exp (-nA^-^/^^^y 



148) 



n=l 



where 



S, 



Si, if ReSi > 0, 
-Si, if ReSi < 0. 



(149) 



Using ( |142|) and (|148|) we see that the conditions (|1 19|) can be written in the form 



-. oo 

= -{z^ + eKA'/^wir/' + y2 



A 



in) 



n=l 



J I 



n=l 



s + nA^-^/^^Si 



(150) 



where Al and M*^"^ are collections of coefficients and 

MW = (z, + ei,Ai/^l;J^-i-^. 



(151) 



The apparent singularity in (|150|) at s = is of course not present provided Mj^^ is 
invertible. M^^-* is of the Vandermonde form and therefore 



detM(°) = (e^Ai/^) 



(K-l)iK-2)/2 



p>q 



Wp - Wq) 



(152) 
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Since none of the roots are degenerate M*^°^ is invertible. Therefore the first non- 
analyticity occurs at s ~ —^^^^/"^^ and hence the mass gap exponent is 



-=l-^ (153) 



6 Discussion 



We have found that v = 1 — 3/2i^ and rj = 2 + 2/{2K — 3) and therefore, since 
7 = the Fisher scahng relation (|13]) is indeed satisfied. The results show that 

the functions which are initially undetermined in the peeling calculation do not in 
fact change the conclusions one would draw simply by ignoring their contribution 



in (p6|). The value of r] agrees with that obtained by slicing [pT^ . 

To compare v it is necessary to examine the continuum limit for the perturbation 
( |132| ) away from the multi-critical point. The multi-critical models have K inde- 
pendent coupling constants so there are K independent directions along which the 
multi-critical point can be approached; if chosen appropriately these directions cor- 
respond in the continuum limit to operators of (length) dimensions {1,2,3, . . .K} 
|T^. The procedure for taking the continuum limit is explained in By definition 



boundary loops with / legs in the dual graph (ie which are / links long) have contin- 
uum length L = la where a is the length of one link. This implies that the generating 
function variable x conjugate to / must be related to a continuum quantity X by 
X = XcC'""^ where Xc is the radius of convergence of the disk amplitude ( |138| ). We 
can construct a non-trivial continuum disk amplitude Ac{A,X) provided that 

A ~ Aa^ (154) 

where A is some continuum coupling; we obtain 

A{A, x) = Ao{A, x) - a^-^/^Ac{A, X) + . . . (155) 

where Aq is the analytic (non-universal) part of the disk amplitude. Note that A 
is not the cosmo logical constant; in the continuum theory the volume V must have 
dimension (length)^ and therefore the cosmological constant A must correspond to 
whichever direction in the space of lattice couplings leads in the continuum limit to 
an operator of dimension 2. It is straightforward to check that this is accomplished 
by the modified couplings 



^ ^ {K-p-2)\2p\ 



9^ = 9^. p = K,K-l. (156) 
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The weight for polygons of 2p sides is modified by a p-dependent factor and so the 
lattice quantity whose continuum limit is the volume is a complicated object with 
polygons of different number of sides weighted in different ways; the volume is not 
the number of polygons. 

Using the scaling ( |154D for A in the large t behaviour of the two-point function 



we find that 

exp(-A'^t) = exp(-A^(to^^)) = exp(-A^T) (157) 

and hence we will get a consistent non-trivial scaling limit provided the continuum 
string time scales as T ~ ta^^ . This also agrees with (where what we call the 
string time is called the geodesic distance). 

The structure of these multi-critical surfaces seems slightly bizarre. Recall that 
the coupling A is conjugate to the number of polygons which therefore behaves 
roughly like A~^ for typical surfaces in the ensemble. On the other hand from (|157| ) 
we have that the typical t must be roughly A~^ and therefore 

(^polygons) ~ (t) . (158) 

When K gets large v 1 and so these surfaces have an almost linear structure when 
viewed in terms of polygons. However, as we discussed in section 1, the polygons can 
be very large so that, for example, the number of links traversed in a typical cycle 
can also be very large and it does not follow that the surfaces in the continuum limit 
are at all one-dimensional. It does however cast doubt on the idea that it makes 
sense to identify the string time T with a continuum geodesic distance. 

Another manifestation of this can be seen by considering a perturbation in terms 
of the cosmological constant A. Then we expect that the two point function behaves 
as 

exp(-A^T) (159) 

(it is a straightforward calculation along the lines of section 5 but using the couplings 
( |156| ) to check this). The average volume of the system will be {V) ~ A~^; on the 
other hand from (|159|) the string time extent must be of order Tg ~ A^"^ [| and 

2 

hence (y) ~ Ts'^" . Interpreting T as a geodesic distance leads to the conclusion that 
the Hausdorff dimension is dn = S- which vanishes as i^' ^ oo and is less than 
two for all K > 2. However we can look at this another way; in the probability 
distribution for the length L of exit loops at time T given a point-like entrance loop 
at T = was calculated. At large but finite K the distribution shows that there 

'''One might worry that the negative weights appearing in the multi-critical models cause some 
cancellation which alters this conclusion. We have checked explicitly that for K = 4 a.t least this 
does not happen. 
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is typically one macroscopic exit loop of length L ~ T^k^ but that as K ^ oo the 
distribution function becomes a delta-function. Using this result we can relate the 
volume to the typical boundary length as K — > cx) and find 

V L'^ (160) 

which is typical of a smooth flat surface of HausdorfT dimension 2. This is the 
behaviour we expect because K ^ oo corresponds to central charge c —oo 
where semi-classical properties are recovered; the volume and loop length have highly 
anomalous behaviour relative to the string time but exactly what we expect relative 
to each other. Note that our discussion above implies that the macroscopically large 
boundary loops are often made from a finite number of polygons with diverging 
number of sides. 

The multi-critical models are non-unitary and have negative central charge, and 
so it is plausible that in the continuum theory the Hausdorff dimension is given by 
(1^ which works well for the c = —2 model. As c ^ — oo gives dn = 2 and so it 
is clear that the geodesic distance implicit in the derivation of this formula (which is 
done directly in the continuum using scaling arguments and Liouville theory) cannot 
be equivalent to the continuum limit of the string time, T. Since we do not know 
what the relationship between these distance measures really is, it is not clear what 
calculation in the discretized formulation would be needed to check (|^) for general 
K] in any case it must be highly non-trivial because (|^) gives irrational values for 
du when K > 2 whereas discretized calculations are almost sure to give rational 
values for exponents. 

We acknowledge useful conversations with Thordur Jonsson. This work was 
supported in part by PPARC Grant GR/L56565. 
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